The scheme for the numerical solution of the incompressible Navier-Stokes equations coupled with the equation for temperature through the temperature dependent viscosity and thermal conductivity coefficients is presented. It is applied, together with the spectral element method, to the 2D calculations of flow around heated cylinder. High order polynomial approximation is combined with the decomposition of whole computational domain to only a few elements. Resulting data are compared with the experimental...